hc/Ae oscillations in a model for (100)/(110) SQUIDs of d-wave superconductors 
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We use a model of hard-core bosons to describe a SQUID built with two crystals of d^2_y2- 
superconductors with orientations (100) and (110). Across the two faceted (100)/(110) interfaces, 
the structure of the superconducting order parameter leads to an alternating sign of the local 
Josephson coupling, and the possibility of quartet formation. Using a mapping of the boson model 
to an XXZ model, we calculate numerically the energy of the system as a function of the applied 
magnetic flux, finding signals of hc/Ae oscillations in a certain region of parameters. This region 
has a large overlap to that at which binding of bosons exists. 
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I. INTRODUCTION 

The interference between two superconducting conden- 
sates results in a wide range of interesting phenomena. 
In particular, it has been suggested that in the cuprates, 
the Josephson effect provides a direct probe of the angu- 
lar dependence of the superconducting order parameter 
[1], and several type of phase sensitive experiments have 
been performed, which confirm the (i2,2_y2-symmetry of 
the superconducting state in high- Tc cuprates [2] . 

If two crystals of cuprate superconductors are con- 
nected at an interface, which is perpendicular to the (100) 
or (010) direction in both crystals, one can have a conven- 
tional Josephson junction (0-junction) or one with a sign 
reversal of the Josephson coupUng (so called 7r-junction) 
depending on the mutual orientations of the (ij,2_y2-wave 
superconducting order parameter [1] and [3]. In prin- 
ciple (neglecting multiple Andreev reflections [4]), for a 
perfectly flat (100)/(110) interface between two crystals 
of cuprate superconductors, the Cu02 lattices meet at 
45°, such that a lobe (antinodal direction) of the d^'^-y^- 
order parameter of one superconductor points towards 
a nodal direction of the other, and therefore the Joseph- 
son current vanishes by symmetry. However, microscopic 
roughness allows for local Josephson tunneling across the 
interface facets, with the sign of the coupling depending 
on the orientation of each facet [5] . This leads to a variety 
of interesting effects like spontaneous supercurrent loops 
[5], locally time-reversal symmetry breaking phases [6], 
[7] and [8], or anomalous field dependencies of the criti- 
cal current density [9]. 

A particularly interesting experimental observation in 
superconducting quantum interference devices (SQUIDs) 
with (100)/(110) interfaces is the hc/Ae periodicity of the 
critical current with applied magnetic flirx [10]. This cor- 
responds to periodicity of half a flux quantum. Since the 
natural explanation for the usual periodicity of one flux 
quantum $o = hc/2e is that the object that tunnels (a 
Cooper pair) has charge 2e, a possible explanation of 
the periodicity of hc/Ae is that electrons tunnel in quar- 
tets, with total charge 4e. Although at first sight this 
possibility seems exotic, it has been proposed before in 



nuclear physics, where a formation of a four-particle con- 
densate was predicted at low density [11]. Furthermore, 
it has been proposed that pairing-quartetting competi- 
tion is expected to be a general feature of interacting 
fermion systems [11]. Motivated by the dominance of 
the second harmonics (periodicity of hc/Ae) of the de- 
pendence of the current with flux, Hlubina et al. [12] 
studied a model for fluctuations of the phase in an array 
of Josephson junctions. They conclude that the cuprates 
are close to an exotic phase with quartet condensation. 
Furthermore, while in principle one might think that a 
more conventional explanation of half flux periodicity is 
the vanishing of the first harmonic in a phcnomcnologi- 
cal Hamiltonian for Josephson junctions [10], the work of 
Hlubina et al. suggests that this is related with quartet 
formation. Thus, the microscopic origin of the domi- 
nance of the second harmonic seems to be the quartet 
condensation. 

Recently, the alternating sequence of superconducting 
0- or TT-junctions (corresponding to a faceted (100)/(110) 
interface) was modeled by a bosonic lattice Hamiltonian 
with a staggered sign for the hopping amplitude across 
the facets [13]. Each boson represents a Cooper pair and 
the boson hopping through the facet is the Josephson 
coupling energy. As a consequence of partial frustration 
of the kinetic energy, the tendency towards boson pair 
formation in the presence of a weak attractive interaction 
is strongly enhanced and in some region of parameters 
quartet formation is favored at the interface [13]. As 
a consequence hc/Ae oscillations are expected in closed 
loop with (100)/(110) interfaces, such as the SQUIDs of 
the experiments of Schneider et al. [10]. The source 
of attraction between pairs might be of the same origin 
as the magnetic mechanism believed to be the source of 
binding in the cuprates. In fact simple arguments used 
to justify attraction of holes in an antiferromagnet, or a 
system with short range antiferromagnetic order, can be 
easily extended to more particles [13]. However, usually 
the loss of kinetic energy does not allow binding of more 
than two particles. 

Since Josephson-junction-type Hamiltonians can be 
derived from a boson lattice model when fluctuations of 
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the amplitude of the superconducting order parameter 
are integrated out [14] and [15], as stated above, one 
might expect that the approaches of Refs. [12] and [13] 
are related and that the appearance of higher harmon- 
ics in the Josephson current is directly connected with 
quartet formation. 

In this work, we calculate the energy as a function of 
the flux £'($) for a simple model describing the motion 
of Cooper pairs in a SQUID with two (100)/(110) inter- 
faces, as in the experimental situation [10] (see Fig. 1). 
Since the current through the ring is proportional 
to dE/dip for a given geometry, the oscillations in E{^) 
are directly related to those of the current. While some 
ideas were reported before [13], the calculation of pair- 
ing was limited to the dilute case of very few bosons in 
the system, and _E((f>) was not calculated. In Section 2 
we explain the model used and its mapping to a spin 
1/2 model. The results of the numerical solution of this 
model are presented in Section 3. Section 4 contains a 
short summary and discussion. 




FIG. 1: Sketch of a SQUID with two faceted (100)/(110) in- 
terfaces and the simple Hamiltonian Eq. (1) used to describe 
it. 



II. THE MODEL 

We describe the Cooper pairs as bosons moving on a 
lattice [13]. At the interfaces, the boson hopping across 
a facet represents the Josephson coupling energy of the 
quantum phase Hamiltonian. 

The Hamiltonian is [13] (see Fig. 1) 

H = H1+H2 + H3 

Hi = Y,(-''^i-'^yaii+l"■o,■^~t'"■i^+2'^c,^+H.C.), 
ai 



ai 

H3 = -^i5](a^^.a2,e^^(-'^''''+i^.c.) (1) 

j 

where al^- creates a boson (Cooper pair) at site i of the 
interface a— 1,2. The first term Hi describes the kinetic 
energy of the bosons within each interface. The hopping 
with amplitude t describes Josephson tunneling across 
the (100)7(110) interface and has an alternating sign due 
to the roughness of the surface and the different change 
of phase of the superconducting order parameter. The 
term in t' describes motion along the interface without 
crossing it. The second term H2 represents an on-site 
U and a nearest-neighbor V interaction. The last term 
H3 describes in the simplest possible way the motion of 
Cooper pairs from one interface to the other. The phase 
(p is related to the magnetic flux threading the ring by 
(p — 27r$/$o, where $0 = hc/2e is the flux quantum. 
For simplicity we take U +00. This implies hard- 
core bosons and allows to map the model into a spin- 1/2 
model as described below. The nearest-neighbor inter- 
action is assumed attractive {V < 0) unless otherwise 
stated. Its origin can be the same as the magnetic pair- 
ing interaction in the cuprates [13]. From results of the 
equivalent spin-1/2 XX Z model for one interface [16], 
we expect that inclusion of second nearest-neighbor at- 
traction does not change essentially the physics. The 
values of the parameters were estimated as V ^ —250 
K, t' > t ^ 100 K [13]. These parameters favor quartet 
formation at one (10G)/(11G) interface, but not in con- 
ventional interfaces with 0-junctions only [13]. 

The alternating sign in the first term of Eq. (1) can be 
eliminated through an adequate change of phases of the 
boson operators, 

— — — — 

which transforms the first term into 

Hi = Y.[-tbl+ib^, + t'bl^^b^, + H.c] (3) 

ai 

while the other terms retain the same form with the re- 
placement aai — > bai- 

Using a boson-spin transformation 

St^ - S- = (-1)'6„., 5^, = blb^, - i (4) 

the model is mapped into the following spin-1/2 model 

Hs = ''^^[JliSai^ai+l + + ^^ai^ai+l) 

ai 

+ J2iSaiSai+2 + ^ai^ai+2)] 

+ J3E(^''^^^'^'''^5^2,- + ^-^-)/2 (5) 
i 
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with the spin exchange couphng constants Ji = 2t, J2 = 
2t', J3 — 2t, and the anisotropy parameter A = V/2t < 0. 
For J3 = 0, the model describes two uncoupled XXZ 
chains representing the interfaces. If periodic or antipcri- 
odic boundary conditions arc used for each chain, the 
Hamiltonian Eq. (5) can be thought as describing two 
XXZ rings, one on top of the other with a complex spin- 
flip interaction J3 between them, the phase of which al- 
ternates sign between odd and even sites of the rings. The 
spin- 1/2 XXZ chain has been studied before in the con- 
text of metamagnetic transitions and results for magnon 
binding (bosons in the original language) and phase sep- 
aration were obtained in specific parameter regimes [16]. 
The region of parameters for which binding of two bosons 
occur, and clustering in more than two bosons or phase 
separation takes place has been studied in more detail in 
Ref. [13]. Taking t = 1 as the unit of energy, the re- 
gion most favorable for pairing of bosons can be defined 
roughly as t' > 1 and 1 < -V < 2t' + 2. In the fol- 
lowing section, we report our search for signals of hc/Ae 
oscillations in E{^) in parameter space and for a finite 
concentration of bosons. 

For this model, the current in the loop is given by [17] 



jm = 



2e d{H) _ 2eNs dE{^) 



h dip 



dtp 



(6) 



where Ns — 2L is the number of sites in the system, and 
L is the number of sites in one ring. 



III. RESULTS 

We have studied by numerical diagonalization the 
equivalent spin Hamiltonian Eq. (5) for two rings of L 
sites each and a number N of up spins (bosons in the 
original language) in a background of down spins. The 
total spin projection Sz = N ~ L/2 is a conserved quan- 
tity in Eq. (5). In general, we have taken 8 < L < 12. 
For small N some calculations were done also for i = 14 
(28 sites counting both rings). We have chosen for each 
ring the boundary conditions (BC) periodic or antiperi- 
odic, which lead to the minimum ground state energy per 
site as a function of flux E{^). In all studied cases, these 
BC did not change as the flux $ was varied. Actually in 
general, and particularly in the region of binding, the BC 
which lead to the minimum energy are twisted [16] and 
[19]. However, the deviation from periodic or antiperi- 
odic BC is small for t' > t and we believe that our choice 
does not affect the results. 

The energy per site i?($) is periodic in one flux quan- 
tum: £'($) = £'($ + $0)- We have investigated the 
tendency towards a periodicity of half a flux quantum in 
E{^). This is reflected in the presence of two relative 
minima in E{^) for <i> = and $ = <&o/2. For N ^ 2 
and = 4, we have explored the region of parameters 
--3 < t' /t < 3, —8 < V/t < 0, examining in more detail 
the region of binding determined previously [13] (roughly 
t'/t > 1 and 1 < -V/t < 2t' + 2,t > 0). For only two 



bosons in the system {N = 2), the region of binding 
has been determined analytically in the thermodynamic 
limit [13] and [16]. Rather surprisingly, we do not find 
two minima for TV = 2 in all the explored region of pa- 
rameters. Instead, for = 4, we do find clear signals of a 
periodicity in <I>o/2, and only inside the region of binding. 
An example is shown in Fig. 2 for t'/t > 1 and L = 8. 
Similar results are obtained for larger system sizes and 
moderately larger values of t'/t. However, keeping A^ = 4 
and increasing L up to L = 14, the relative minimum of 
i?($) at $ = becomes less pronounced, suggesting that 
it disappears and there is no periodicity of E{^) in ^o/2 
in the thermodynamic limit for dilute systems N/L — > 0. 
In any case, the dilute limit is not expected to be rele- 
vant to the experimental situation. To be able to analyze 
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FIG. 2: Ground state energy per site as a function of the 
applied flux for L = 8.N = 4, t'/t = 1.5, V/t = -2 and t±/t = 
0.1. 

how E{^) evolves with system size at a finite fixed con- 
centration of bosons, we calculated £'($) for increasing 
L, keeping N — L and the other parameters fixed. Note 
that due to the symmetry of the spin Hamiltonian Eq. 
(5) under a change of sign Sz, the problems with A^ or 
2L — N bosons are equivalent. Therefore, the case with 
N — L is that of the greatest possible total kinetic energy 
in absence of the interaction V. In Fig. 3 we show how 
E{^) changes as L increases from 8 to 12 (in even steps 
to avoid frustration) and other parameters as in Fig. 2. 
The absolute minimum is at $ = and there is a relative 
minimum at = $o/2. This minimum is rather shallow 
for L = 8, but becomes more pronounced as L increases. 
This is consistent with a periodicity of E{^) in half a fiux 
quantum in the thermodynamic limit. For a reasonable 
value of the Josephson coupling energy t ^ 10~^ eV, one 
has Aet/h ~ lOfiA. Since the number of facets in a real 
system is of the order of L = 10^ - 10"* [18] , from Eq. (6), 
the current extrapolated to the thermodynamic limit is 
j ~ 1 — lOfiA. This in rough agreement with the exper- 
imental order of magnitude j ~ 10 — 100^^ [10]. Since 
roughly j ~ i^(see below) the agreement improves if t± 
is increased. 

The evolution of E(^) with size for larger t' is shown 
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FIG. 3: Ground state energy per site as a function of the 
applied flux for t' /t = l,V/t = -2,t±/t = 0.1, iV = L and 
different system sizes L. 



in Fig. 4. In the spin language, these parameters cor- 
respond to a next-nearest-neighbor antiferromagnetic in- 
teraction in the x and y directions J2 two times larger as 
the corresponding one for nearest neighbors Ji. In this 
situation, in the classical case, an antiferromagnetic order 
between next-nearest neighbors is favored and to avoid 
frustration of this order in a ring, the number of sites L 
should be multiple of 4. As a consequence, although this 
frustration is only partial in the quantum case, we believe 
that the results for iV = L = 10 are not reliable for t' > t. 
This is the case of the middle panel in Fig. 4. Comparing 
the other two cases represented in Fig. 4, although none 
of them shows two well defined minima, the results are 
not inconsistent with a double periodicity in the thermo- 
dynamic limit because has a rather sharp maximum 
for Af = L = 8at$ = 0.5$ that becomes rather flat for 
N = L ~ 12, suggesting the development of a minimum 
as L is further increased. Note that for t±/t — 0.2, a 



relative minimum in E{^) at $ = 0.5$o already exists 
for TV = L = 12 (dashed hue in Fig. 4). We have also 
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FIG. 4: Same as Fig. 3 for t' /t = 2. The dashed line at the 
bottom is the result for t±/t — 0.2 multiplied by a factor 0.2. 

analyzed the effect of varying hopping between interfaces 
tj_. As t± increases, the tendency towards an additional 
periodicity increases slightly up t±/t ~ 0.2 and for larger 
tj_ , it is weakened. This is expected, since the additional 
kinetic energy tends to break the bound pairs. For small 
t±, from perturbation theory, the structure of E{^) is 
dominated by several terms of order t']_ . This is in rough 
agreement with our numerical results for t±/t ~ 0.2, al- 
though the increase in the amplitude of E{^) with t± 
seems faster than quadratic (see for example the bottom 
of Fig. 4). 

Most of the above-mentioned results correspond to the 
half filled case N = L. We have also looked for signals 
of an extra periodicity in i?(<&) in the quarter filled case 
N — L/2, comparing the flux dependence of the energy 
for L = 8 and L = 12. The results suggest that no double 
periodicity is present in the thermodynamic limit. 
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Finally, we also study the shape of for repulsive 

interaction > in the half filled case. Surprisingly, 
we find a tendency for double periodicity for positive t' 
(see Fig. 5). We have studied the binding energy of an 
isolated ring Ab = L[EiiN + 2) + Ei{N)~2Ei{N + l)] for 
different parameters, where Ei{N) is the energy per site 
of one ring for N = L/2 particles. The region of binding 
coincides roughly with that obtained previously for N — 
[13] and [16] and > for positive V. The absence of 
a negative binding energy for positive V also persist when 
t± is included. The factor of 10 increase in the magnitude 
of the oscillations as L is increased form 8 to 12 points 
to particularly large finite-size effects. In all cases we 




FIG. 5: Ground state energy per site as a function of the 
applied flux for t' /t = 1.5, V/t = -2,t±/t = 0.2 and different 
values of N = L: full line L — 12, dotted line L = 8, dashed 
line 10 times L = 8. 

do not find signs of extra periodicity for negative t' . We 
remind the reader that negative t' corresponds to the case 
of ordinary junctions in the facets (all 0-junctions), for 
which no extra periodicity was experimentally observed. 

IV. SUMMARY AND DISCUSSION 

We have calculated numerically the energy per site as 
a function of applied magnetic flux E{^) in a simplified 
model for hard-core bosons (representing Cooper pairs) 
for a SQUID containing two (100)/(110) interfaces. In 



spite of the limitations of the size of the systems studied, 
the finite-size scaling suggests hc/Ae periodicity for pa- 
rameters for which binding of bosons at one interface take 
place [13]. This fact and the magnitude of the resulting 
current is consistent with the experiments of Schneider 
et al. [10]. This periodicity is absent for usual interfaces 
it' < in the model). 

We also find signals of hc/Ae oscillations in E{^) for 
repulsive interactions for which no binding exists. The 
reason of the extra periodicity in this case is unclear, and 
needs further study. 

For simplicity the motion of Cooper pairs inside one 
crystal from one interface to the other has been repre- 
sented by only one hopping parameter t±. This is a 
rather crude description which allowed us to include more 
sites at the interfaces, keeping the total amount of sites 
below 32 due to computer limitations. In any case, we 
expect that the essential physics is retained, at least for 
attractive interaction between Cooper pairs. The phys- 
ical picture can be the following: quartets (bound pairs 
of bosons) are formed at the interfaces as a consequence 
of the competition between attraction and the reduced 
kinetic energy there [13]. The energy scale of the mo- 
tion perpendicular to the interface {t± in our model) is 
responsible for the coherence and the observed flux de- 
pendence. However, if it is too large, it tends to unbind 
the quartets. If this image is correct, and quartets live 
only near the interface, the observation of hc/Ae oscil- 
lations should depend on the size of the system. They 
should be weaker or disappear if the coherence length in 
the direction normal to the interface is much smaller than 
the size of the system. 

A more realistic model for the description of the mo- 
tion of Cooper pairs between both interfaces should in- 
clude one or more layers of intermediate sites between 
the interfaces. 
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